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Abstract. We construct Backlund transformations (BTs) for the Kirchhoff top by ta- 
,_^ king advantage of the common algebraic Poisson structure between this system and the 

T-H sl{2) trigonometric Gaudin model. Our BTs are integrable maps providing an exact time- 

^D discretization of the system, inasmuch as they preserve both its Poisson structure and its 

^^ invariants. Moreover, in some special cases we are able to show that these maps can be 

^H explicitly integrated in terms of the initial conditions and of the "iteration time" n. En- 

. ^ couraged by these partial results we make the conjecture that the maps are interpolated by 

a specific one-parameter family of hamiltonian flows, and present the corresponding solu- 

^ ' tion. We enclose a few pictures where the orbits of the continuous and of the discrete flow 

I— I are depicted. 
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>. 1 Introduction 

^^ The Kirchhoff top is an integrable case of the Kirchhoff equations [1] describing the motion 

pg of a solid in an infinite incompressible fluid. In general the total kinetic energy of the system 

1^ solid + fluid is given by a quadratic expression both in the translational velocity u of the 

(^ rigid body relative to a fixed frame and in its angular velocity u [2]. If the solid has three 

"O perpendicular planes of symmetry and is one of revolution too, say around the z axis, or is 

. . a right prism whose section is any regular polygon, then the total kinetic energy reduces to the 

^^ simple diagonal form [3]: 

X 1 1 

"y T = -{A.iuj + uD + A^ul) + -{Bi{ujJ+ujI) + B3uI), (1.1) 

where the quantities Ai, A3, Bi, B3 are constants depending on the particular shape of the 
solid. The total impulse p and angular momentum J of the system, i.e. the sum of the impulse 
and angular momentum of the solid and those applied by the solid to the boundary of the fluid 
in contact with it, are given by [2]: 

_ dT J _ dT 

Pi — T^ ) 'Ji — 7^ ■ 

OUi OUJi 

By an Hamiltonian point of view, impulse and angular momentum must obey the Lie-Pois- 
son e(3) algebra given by the following Poisson brackets: 

{Ji, Jj} = CijkJk, {Ji,Pj} = eijkPk, {Pi,Pj} = 0. (1.2) 



*This paper is a contribution to the Proceedings of the Conference "Integrable Systems and Geomet- 
ry" (August 12-17, 2010, Pondicherry University, Puducherry, India). The full collection is available at 
http://www.emis.de/journals/SIGMA/ISG2010.html 
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where i, j, k belong to the set {1, 2, 3}. These brackets have two Casimirs: 

3 3 

Y^mJi = Ci, Y.pj = 2C2. (1.3) 

1=1 i=l 

Rewriting the kinetic energy (1.1) in terms of the piS and Jj's, one has two commuting integrals 
of motion for the Kirchhoff top: 



1 (pI + pI , pIA 1 (^1+4 , 4 



The flow with respect to the Hamiltonian T is given by the expressions 

p = {r,p}, J = {r,j}. (1.4) 

2 The Kirchhoff top by a contraction 
of trigonometric Gaudin model 

In this section we show how to obtain the Lax matrix for the Kirchhoff top, in all the cases when 
the relation B^ = A^ — A^ holds, by a procedure of pole-coalescence on the Lax matrix of 
the two-site trigonometric Gaudin model [4]. The main results are derived in [5, 6]. To this aim, 
let us briefly review some relevant features of the trigonometric Gaudin model. In the two-spin 
case the Lax matrix reads: 

( Ag{\) Bg{X) \ ,_. 

2 2 — 2 + 

AaW=Yj-^(X-X,)Jl, Bc(A)=E-^. £■„(>) = ^-^^^ (2.2) 

In (2.1) and (2.2) A G C is the spectral parameter, Aj are the arbitrary parameters of the Gaudin 
model, while (s^, s^, s^), j = 1, . . . , 2, are the spin variables of the system obeying to ®^s/(2) 
algebra, i.e. 

I^i' 4} = T^^ifc^fc > 14' «fc } = -2^'^ifc4- (2.3) 

In terms of the r-matrix formalism, the Lax matrix (2.1) satisfies the linear r-matrix Poisson 
algebra: 

{LG(A),LG(/i)} = [ri(A-^),LG(A)»/ + I»LG(/i)], 
where rj(A) stands for the trigonometric r matrix [7]: 



niX) 



sin(A) 



/ cos(A) \ 

10 

10 

V cos(A) / 



The determinant of the Lax matrix (2.1) is a generating function of the integrals of motion. In 
fact we can write: 

A^ifT (\\\- ^^Q I '^^g I iJ'Gsin(Ai-A2) 2 

sm(A — Ai)^ sm(A — A2) sm(A — Aij sm(A — A2) 
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where Cig and C2G are the Casimirs of the algebra (2.3) given by Cjg = (sf)^ + s^s- , while 
the two involutive integrals of motion Hq and Hq are: 

o- _ 2C0S(A1 - X2)slsl + Sfs^ + S^st o- _ ,3 , ,3 ^ t3 j rr rr \ - H 

J^G — r-7T r^; , J^O — Si + S2 — Jg, i-nG,-noj— U. 

sm(Ai — A2J 

To get the Kirchhoff top we perform the pole-coalescence by introducing the contraction pa- 
rameter e and take in the Lax matrix (2.1) Ai — t- eAi and A2 -^ eX2- The Lax matrix for the 
Kirchhoff top is recovered by setting: (the notation is v^ = v] it ivf, Vj = {vl,vf,vf) for any 
vector set Vj): 

J = si + S2, p = e(AiSi + A2S2) (2.4) 

and letting e — )• in (2.1) after this identification. By using (2.3), it is readily seen that the 
variables J and p (2.4), obey the Lie-Poisson algebra e(3) (1.2). Finally, the Lax matrix for the 
Kirchhoff top reads: 

^3 



, A{X) B{\) 
^y^) - 1 c{X) -A{X) 



I /xN X-:! V J COt(A)» \ 

/ cot A J^ + — ^ —1V\+ ■ (\\ 1 

sm(A)"' sm(A) sm(A) 

J+ cot(A)p+ , ,,. -,0 P^ X 

V -TYT+ ■ M ^ -cot A -73 + — ^ j 

\ sm(A) sm(A) sm(A)^ / 



(2.5) 



Again, its determinant is the generating function of the integrals of motions. Indeed we have: 

-det(L(A)) = ^^ + 2//ocot(A)2 + 2C2^^ + 2Ci^^, (2.6) 

sm(Aj^ sm(Aj^ sm(Aj^ 

where Ci and C2 are the Casimirs (1.3), while Hq and Hi are the two commuting integrals 
given by: 

Hi = \{Jl + Jl+pfl, 2Ho = Ji, {Hi,Ho} = 0. (2.7) 

In all cases where B^ = A^ — A^ , the total kinetic energy (1.1) can be rewritten in terms of 
the quantities (1.3), (2.7): 

,„ G2 tin ll^ , , 

3 Backlund transformations 

In this section we construct a two parameter family of Backlund Transformations defining sym- 
plectic, integrable and explicit maps that, as we will see, provide an exact time-discretisation 
of our model. The approach follows that given for example in [8] and take advantage of the 
results derived in [9] where the Backlund transformations (BT) for the A^-site trigonometric 
Gaudin magnet have been constructed. In fact, since the r-matrix structure survives the pole- 
coalescence and contraction procedures, the ansatze for the dressing matrix -D(A) linking, by 
a similarity transformation, the old Lax matrix L{\) to the new Lax matrix L{\) are the same 
as for the trigonometric Gaudin. Thus, according to the procedure followed in [9], we write: 

L{\)D{\) = D{X)L{X), (3.1) 

where L has the same A dependence as in (2.5) but is written in terms of the updated variables 
(J^ J=^,p^p=^). The matrix D{X) reads [9] 

r^/>N ^ / sin(A- Ao -/u) + P(5cos(A- Ao) Pcos(//) \ , . 

^ ' I Q sin(2/u) - PQ2 ^^^g^^) g^^^^^ _ Aq + ^) - PQ cos(A - Aq) T ^ ' 
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In (3.2) Ao and fi are arbitrary constants and P and Q are, up to now, indeterminate dynamical 
variables. 

We remark that in the fundamental paper by V. Kuznetsov and P. Vanhaecke [10], an exten- 
sive study of BTs in the 2x2 case has been performed. In fact a reformulation of our similarity 
transformation (3.1) having a polynomial dependence on the spectral parameter can be derived 
starting from their results. 

Our aim is now to find an expression for P and Q in terms of only one set of dynamical 
variables, say the old ones, so that (3.1) yields the explicit map between the two sets of variables. 
To achieve this goal, we use the so-called spectrality property (see for example [8]). 

Note that the determinant of D{X) is proportional to sin(A — Xq — ^) sin(A — Aq + /u), so, 
modulo 27r, it has two zeros, A+ = Xq + h and A_ = Aq — /u. D(X±) are clearly rank one matrices, 
having one dimensional kernels, say, \K±). The key point is that these kernels are eigenvectors 
of the Lax matrix. Indeed from (3.1) it follows: 

L{X±)\K±)=j±\K±), (3.3) 

where the two eigenvalues are given by: 



^l = A\X) + B{X)CiXy 



A=A± 



and ^(A), B{X) and C(A) are defined in (2.5). The equation (3.3) gives the relations be- 
tween P, Q and the old dynamical variables. In fact, the two kernels are given by: 

'^+^=(-q)' '^-^ = (2sin(^)-PQ 

and then readily follow the expressions for Q and P: 

n n(x \ ^ Q(A+)-Q(A-) A{X±)ti{X±) 
^ = ^^^+)' P = 2 sin(/.) ' ^(^^) = B(X^) • 

Taking the residue of (3.1) at the pole in A = and its value at A = | we obtain the explicit 
maps as below: 

p" = > . .. ^ ■ ., X (a+P~ - -P^cos(^)^p+ + 2Pcos(^)a+p^), 
Asm(A4.)sm(A_) 

p+ = ........ A a-P'^ - Q'^cos{nfc^p' - 2Qcos(/i)ca_p^), 

Asm(A+)sm(A_) 

\ 

P^ = A ■ f\ \ ■ f\ j 2a+a,/-Pcos(/x)a-p+ + Qcos(//)ca+p~) -p^ 
Asm(A+)sm(A_) ^ 

J- = -^ nr-^i^+'^' - P'^os{fifj+ - 2Pcos(m)6+/), 

Acos(A+j cos(A_) 

J^ = -^ 7V-^ T^^(blJ^-Q^cos{fifc^J- +2Qcos{fi)cb-p^), 

Acos(A+J cos(A_j ' 

J^ = J'-^, (3.4) 

where 

a± = sin(A-i-) =F PQ cos{Xo), b± = cos(A-i-) it PQcos(Ao), 
A = 1 - 2PQ sin(;u) + P^Q^ c = 2 sin(/u) - PQ. 

Thus the maps depend on two Backlund parameters, Aq and ^ (or A+ and A_): in the next 
section we will show that, provided Aq G M and n G iM., this two-point transformation is actually 
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a time discretization of a one parameter family of continuous flows having the same integrals of 
motion (1.3), (2.7) as the continuous dynamical system ruled by the physical Hamiltonian (2.8). 
With the above constraints on the parameters, the BTs become physical, mapping real variables 
into real variables. Furthermore these transformations are symplectic. In fact, as the r-matrix 
structure underlying the Kirchhoff top is the same as that of the ancestor trigonometric Gaudin 
magnet, the simplecticity of the transformations (3.4) is guaranteed along the lines given in [9]. 
Next, we will formulate the conjecture that, provided Aq G M and fi S iM, this two-point 
transformation is not only a time discretization of a one parameter family of continuous flows 
equipped with the same integrals of motion (1.3), (2.7), but it has also the same orbits as the 
continuous dynamical system ruled by the physical Hamiltonian (2.8). The above conjecture 
will be verified to hold in a couple of special cases, where the explicit solution of the recurrences 
defined by the maps (3.4) will be derived, and shown to be interpolated by the solution to the 
evolution equations for the continuous Kirchhoff top. On one hand this confirms the Kuznetsov- 
Sklyanin intuition that Backlund transformations can be used as a tool for separation of variables 
(see also [10]), on the other hand, in light of [10], these results appear quite natural if one 
thinks that the Backlund transformations are translations on the invariant tori, translations 
that corresponds to some addition formulas for families of hyperelliptic functions. 

4 Continuum limit and discrete dynamics 

As shown in [9], to ensure "reality" of the maps (3.4), one has to require the Darboux matrix D 
to be a unitary matrix (possibly up to an irrelevant scalar factor); this holds true iff \± are 
mutually complex conjugate, i.e. iff Aq is real and ^ is pure imaginary. So we set: 

A+ = Ao + «-, A_ = Ao-«-. 

In the limit e — )• the relations (3.4) go into the identity map. Indeed e plays the role of time 
step for the one parameter (Aq) discrete dynamics defined by the Backlund transformations. 
By following [9] , in order to identify the continuous limit of this discrete dynamics we take the 
Taylor expansion of the dressing matrix at order e, obtaining: 

^' ^ ' 27(Ao) V C{Xo) -A{Xo)cos{X-Xo) J ^ '' 

where the functions ^(A), B{X) and C(A) are given by (2.5), and 7(A)2 = ^(A)^ + B{X)C{X). 
By inserting this expression in the equation (3.1) we arrive at the Lax pair for the continuous 
flow: 

L(A) = [L(A),M(A,Ao)], (4.1) 

where the "time derivative" is defined as L = lim ^^ . 
The matrix M(A, Aq) takes the explicit form: 



M(A,Ao) = -7>n 



/A(Ao)cot(A-Ao) ^^^°^ ^ 



sin(A — A, 



0) 



\ iil!^ -^(A„)c„t(A-A„,^ 



In Hamiltonian terms, the system (4.1) reads: 

LiM) = {7(Ao),i^^j(A)}, i,j G (1,2), (4.2) 
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entailing that the variables p and J of the continuous flow obey the evolution equations: 

P = {7(Ao),p}, J = {7(Ao),J}. (4.3) 

It is clear that the dynamical system given by (4.3) possesses the integrals (2.7), because of (2.6). 
Moreover we have some evidences, that will be reported in the following, that the continuous 
and the discrete system share the same orbits too. 

First of all we note that the direction of the continuous flow that obtains in the continuum 
limit from the discrete dynamics defined by the Backlund transformations (3.4), and that of the 
Kirchhoff top (1.4) with the kinetic energy T given by (2.8), can be made parallel. In fact the 
shape of the orbits are unchanged if one takes an arbitrary C^ function of the Hamiltonian 7(Ao) 
as a new Hamiltonian in (4.2), since this operation amounts just to a time rescaling (for every 
fixed orbit 7(Ao) is constant). Accordingly, we take as Hamiltonian function "''^2 ? where w 
is, so far, an arbitrary constant. The expression (2.6) allows to write the explicit equations of 
motion for a generic function of the dynamical variables J^(p, J): 

t( T^ / ^7(Ao)^ ^, ^A / H^ , gocos(Ao)^ \ 

I 2 J [ sm(Ao)^ sm(Ao)^ J 

This has to be compared with with the equations of motion for the physical Hamiltonian (2.8): 

-^(P,j) = {|^ + f,-^(p,j: 

The two expressions coincide by identifying: 

'"^B'^~B{ sm(Ao)^= '^^ \ (4.4) 

In other words, the physical flow is the continuum limit of the discretized one. Now we make 
the following 

Conjecture. For any fixed Aq, there exist a re-parametrization of e, e — )• T^^, possibly depending 
by the integrals and the Casimir functions, such that T = e + O(e^), and at all order in T the 
continuous orbits of the physical flow interpolate the discrete orbits defined by the Backlund 
transformations, provided that Ao is chosen according to (4.4). 

This is equivalent to say that, for any fixed Aq, via the above reparameterization, Backlund 
transformations form a one parameter (T) group of transformations, obeying the linear composi- 
tion law BTt^^ oBTtj = BTy^+^j, "o" being the composition. Note also that, if the conjecture is 
true, then, since at first order in e (and therefore in T) the flow is ruled by the Hamiltonian 7(Ao), 
one has: 

S,n ^ gnT{7(Ao),.}^ = ^ + nT{^{\,),x} + ^{7(Ao), {7(Ao), x}} + • • • , 

where x" means the n-th iteration of the Backlund transformations with the same parameter T. 
In the Figs. 1 and 2 we report respectively an example of the orbit for the variables (p^(t), 
p'^{t),p^{t)) for the continuous flow ruled by the Hamiltonian (2.8) as given in Appendix A and 
of the corresponding discrete flow obtained by iterating the Backlund transformations. The 
initial conditions are the same and the value of Aq has been chosen so to make the continuous 
limit of the discrete dynamics parallel to the continuous flow of the Kirchhoff top. They overlap 
exactly. In the next section, assuming the conjecture to hold true, we will show a way to find 
the parameter T. There, we will give as well analytic results in two particular cases, where 
the continuous flow is periodic, and not just quasiperiodic. Clearly, these non generic examples 
cannot be invoked to support our conjecture: however, we decided to include them in the paper 
inasmuch as they provide an explicit link between discrete and continuous dynamics. 
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20^--^ 20 



Figure 1. initial conditions: p^O) = 15, ^^(o) = -12.13, ^^(o) = -10, ji(0) = 1, .P{0) = -4, 
J3(0) ^ 3. Moments of inertia: Bi = 1, Bg = scc(0.1)2. 




Figure 2. input parameters: p^{0) = 15, ^^(o) = -12.13, ^^(o) = -IQ, ji(0) = 1, ,P{0) = -4, 
J3(0) = 3, Ao = 0.1, e = 0.1. 



4.1 Integrating the Backlund: special examples 

Let us assume to have a smooth transformation, that we indicate with x = f{x,e), where the 
parameter e plays the role of the time step, such that f{x, 0) = x. By x"' we denote the n-th 
iteration of the map, so that x° = x, x^ = f{x,e), x^ = /(/(x,e),e) and so on. Solving the 
Backlund map amounts to find £" as a function of x, n and e. Now we will show that, under 
given assumptions, there is indeed a positive answer to this question. We will follow a simple 
argument, well known in group theory [11]. 

Suppose to do a transformation from x to 5^ with parameter ei and then another one from x^ 
to x^ with parameter £2- Suppose also that there exist a parameter £3 linking directly x to Sp'. 
As the Backlund are smooth, varying continuously ei or 62 corresponds to a continuous variation 
in £3: the Backlund transformations define €3 as a continuous function of ei and £2, say 63 = 
x{^ii^2)- Now consider infinitesimal transformations: a small change in the parameter e take 
the point 5^ to a near point 5^ + dx^: 



^^ + dx^ 



X 



f{x,e + de). 
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But we can arrive at the same point by starting from x^ and acting on it with a transformation 
near the identity, say with the smah parameter 5e: 

x^ + dx^ = f{x\6e). (4.5) 

The relation between the parameters now reads: 

e + de = x(e,5e)- 
Obviously x(e) 0) = e, so: 



'^^ d6e 



Se = T{e)6e. (4.6) 

<5e=0 



The relation (4.5) tells us that: 

6e = C{x^)de. 



,,.^dn,\6e) 



9S^ 5e=0 

The last expression together with (4.6) gives: 

dy r dX 



This means that there exists a function, say V , such that: 

V{x'^) = V{x)+T =^ y(x") = y(x) + nT. 

Formally we can write this expression as 5" = V^^{V{x) + nT). However, for n = 1 we must 
have x^ = f{x, e(T)), yielding x" = f{x, e{nT)). The continuous flow discretized is simply given 
by x{t) = /(x, e(t)) where x is the initial condition (x(t = 0) = x). 

In the following we will present two particular cases, both corresponding to periodic flows, 
where the Backlund transformations can be explicitly integrated. 

Example 1. Consider the invariant submanifold p = {X, 0, Z), J = (0, Y^ 0). Since now Hq = 0, 
the freedom to have a parameter Aq in (4.2) is just a scaling in time, so we can freely fix it: 
by now we pose Aq = f • With this choice the interpolating Hamiltonian flow discretized by 
the maps (3.4) is given simply hy T-L = yjY"^ + Z^. So, as seen at the beginning of this section, 
in order to have real transformations we pose Ai = ^ + ie and A2 = f — ie. The Backlund 
transformation can be now conveniently written in terms of a single function R oi e, X , Y 
and Z: 

- _ 4i?sinh(6)(i?^ + 1) ^ {R^ + If - 4i?^sinh(6)^ 

(i?2 -1)2+4 COsh(e)2i?2 ^ [R2 _ 1)2 + 4^2 cosh(e)2 ' ^ ' ' 

~_ 4i?cosh(e)(i?2-i) {R^ - If -AR^ coBhjef 

(i?2_i)2 + 4cosh(e)2i?2 (i?2_l)2 + 4^2cosh(e)2 ' ^ ' 

~_{R^ + lf-AR^shi\i{ef^ 4i?sinh(e)(/?2 + i) 



R 



(i?2 _ 1)2 + 4^2 cosh(e)2 (i?2 -1)2 + 4 cosh(e)2/?2 - 

Z - v/(^2 cosh(e)2 - 2C2 sinh(e)2) 
X sinh(e) + Y cosh(e) 



Note that the two constants under square root in the numerator of R are the Hamiltonian 
% = ^/Y"^ + ^2 and the Casimir function C2 = — ^ — . To solve the recurrences (4.8) one has 
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to find e as a function of the parameter T defined in (4.7). To this end we first note that 
so that by the relations (4.7) we have: 



de 



e=0 ^ 



\^=rH^f4.= fn , "' =T. (4.9) 

z 2XY Jo 2XY de J^ y/^^s cosh(e)2 - 2C72 sinh(e)2 

Ah that we have to do now is to perform the integral, invert the result to find e as a function of T, 
then plug the result into (4.8) and replace T by nT: this gives the solution to the Backlund 
recurrences. After some manipulations with the Jacobian elliptic functions we arrive at the 
simple result: 

I sn(r y-^^) 

cn(r, ^) cn(r, ^) 

With this position we can write down the expressions for X", l"" and Z": 

-^ _ 4Rsn{nT)cn{nT){R'^ + 1) {R"^ + l)2cn(nr)2 - 4:R'^sn{nTf 

(i?2 _ l)2cn(nT)2 + 4i?2 + [r2 _ l)2cn(nr)2 + AR^ ' 

~„ _ ARcn{nT){R'^ - 1) {R^ - l)^cn{nT)^ - AR^ 

~ (i22_i)2cn(nr)2 + 4i?2 ~ (i?2 _ i)2cn(nr)2 + 4i?2 ' 
~„ _ (i22 + l)2cn(nT)2 - 4i?2sn(nr)2 4i?sn(nr)cn(nT)(fl2 + i) 

(i?2 _ l)2cn(nr)2 + 4i?2 (i?2 - l)2cn(nr)2 + 4i?2 ' 

^ Zcn(nr) - y/CW^ - 2C72sn(nr)2) 
R = 



Xsn(nr) + Y 

where for brevity we have omitted the elliptic modulus n, ^ in the Jacobian elliptic functions 
"sn" and "en". Note that if we pose in (4.9) 2T = t, that is 

1 sn(i,v^2^) 

cosh(e) = ^^, sinh(e) = — 2' Tj_' 



cnU,^"^ir^ en 



2' n ' ^^^^2' n J 

in (4.8), then we have the general solution of the dynamical system ruled by the interpolating 
Hamiltonian flow T-L = \'Z^ + Y^, that is the value that takes the hamiltonian 7(Ao) (4.3) on 
the invariant submanifold considered in this example for Ao=f . The equations of motion are 
given by UX = -YZ, UY = -XZ, UZ = XY. 

Obviously this general solution coincide with that found by a direct integration of the previous 
equation of motion, i.e. with Z = \/2C2 sn(t + v), X = \/2C2 cn(t + v) and Y = 'Hdn{t + v), 
where the elliptic modulus of this functions is again ^y^ and where v is such that sn(t;) = /^t^ - 

Example 2. In the next example we consider the invariant submanifold p = (x,y, 0), J = 
(0,0,2:). Again, in order to have real transformations we pose Ai = Aq + ie and A2 = Aq — ie 
with Aq and e real. In terms of p = x ±iy, the maps (3.4) become: 

z sin(Ao — ie) — \/ z^ sin(Ao — ie)2 + p^p+ _ 



~+ ^ 

p = - 

2 sin 



zsin(Ao — ie) — \J z'^ ■sva.{\o — leY ^ p p^ 
z sin(Ao + ie) — \J z'^ sin(Ao + ie)2 + p~p"' 

sin(Ao + ie) - ^J z^ sin(Ao + ie)2 + p-p+ + . _ , . 

sin(Ao — ie) — ^ z'^ sin(Ao — ie)2 +p~p+ 
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To find the relation defining the parameter T (4.7), first find the expression of -^ 
2iz cos{\q)p~ 



e=0 



dp 



=0 \/ Z^ SITl{\q)'^ +P+P 

then by using (4.7) one has: 

-y^^^in(AoP^rp+p^ P dp" — T — \/-2^^ sin(Ao)2 + p+p~ 
2 cos(Ao) Jp- izp^ 2 cos(Ao) 

r I cos(Ao + ie) cos(Ao - ie) \ 

Jo y Y^z^ sin(Ao + ie^ + p~^p~ ^ z'^ sin(Ao — ie)^ + p^p~ J 

or more exphcitly: 

, ^^^°^(^o)^ = arcsinh f-^ sin(Ao + .e)^ - arcsinh f-^ sin(Ao - ^ e) 

Vz2sm(Ao)2+p+p~ VVP+P / Vv^^P , 

/ z sin(Ao — ie) — y^^siii(Ao — iep-|-p^+ \ 
y2;sin(Ao + ie) — y^z^ sin(Ao + ie)^ + p~p+ J 

The Backlund transformations (4.10) now take the simple form: 

(2icos(Ao)zT \ _ _, / 2icos(Ao)zr \ , 

= \ P , p = exp = p^, 

y' z'^ sm{\oy +p+p- J y y' z'^ sm{XQ)'^ +p+p~ J 



z = z, 



so that again, as expected, the n-th iteration of the maps {p )", (p"*")"", z" is found by substitut- 
ing T with riT. By posing 2T = t in the previous expressions and returning to the real variables 
X = ^ ^ and y = ^ ~^f , we have the continuous flow: 

( cos(Ao)2; A . / cos(Ao)2; A 

x\t) = X cos 1 + y sm — ^r , 

y Y'2;2sin(Ao)^ +P+P~ y y y^z^ sin(Ao)2 +p+p- y 

/ cos(Ao)z \ . / cos(Ao)z \ 

y\t) = ycos =t — j;sm =t 

\^^ z^'&m.iX^Y ^p^p- y \^^ z'^sm.iX^Y ^p^p- j 

corresponding to the general solution of the continuous system ruled by the value that takes the 
hamiltonian 7(Ao) (4.3) on the invariant submanifold p = (x,2/,0), J = (0, 0, z) considered in 
this example. 

A Integration of the continuous model 

For the sake of completeness, we report here the solution to the evolution equations defining the 
Kirchhoff top. Gustav Kirchhoff in his ^^Vorlesungen iiber mathematische Physik" [1] deals with 
the problem of integrating equations of motion (1.4). For the physical variables (J^, J^, J^) and 
{p^,p'^,p^) they are written in extended form as: 

p\t) = aj\t)p\t) - (3j\t)p\t), p\t) = fiJ\t)p\t) - ap\t)j\t), 

f{t) = p{j\t)p\t) - J\t)p\t)), 

j\t) = (a - P)j\t)j\t) + Pp\t)p^{t), ,P{t) = (/3 - a)j\t)j\t) - f3pHt)pHt), 
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j3(t) = =^ J^ = const = M, 

where for simphcity we have posed a = B^ and /3 = B^ . The new variables suggested by 
Kirchhoff are given by p^ = scos(/), p^ = ssin(/), J^ = o-cos(V' + /), J^ = asin(^jJ + /). In 
terms of these variables the equations of motion can be written as: 

s{t) = —(t{t)p (r) sin(^(r)), a{T) = —s{t)p (r) sin(V'(T)), 

f{T) = s(r)o-(r)sin(V'(r)), 

^{t) = M- p\r) cos(V^(t)) f ^ + ^"j , /(r) = ^/(r) cos(V^(r)) - ^M, 

\s{t) (t{t)J s{t) 13 

where r = (3t, ( ) = -g:^-. By using the constraints given by the Casimirs and the integral Hi, 
i.e. 2Hi = a'^ + (p^)'^, 2C2 = s^ + (p^)'^, sacos{'ip) + Mp^ = Ci, one can readily obtains the 
equation for the evolution of p^ : 

iff = m' - 2C2){{p'f - 2Hi) - {Mp^ - Cif. (A.l) 

At this point Kirchhoff notes that this equation is integrable and that one can obtain by the 
expression oip^ those for the other variables, but soon after he passes to consider the special case 
J^ = J^ = p^ = 0. At our knowledge the first author to integrate this system was G.E. Halphen 
in 1886 [12], also if some authors gives Kirchhoff as reference for the complete integration of the 
equation of motions^. The expression for p'^{t) of Halphen is written in terms of Weierstrass p 
function; it reads as: 



where 

D 2 

Note that the last two arguments of p are not its periods but the elliptic invariants. In order to 
fit with the given initial condition p^{0) it is possible to show that one has to choose the value 
of K according to the relation: 

p{K, $ + 3^^, *3 _ ^$ _ ^73^ ^ \i^±J2 i^J:i^ _ 

By the expression (A. 2) it isn't simple to see that p^{t) is actually bounded. It is however possible 
to make plain this point by writing the solution in terms of the roots of (A.l). Let us clarify 
this statement. The key observation is that if one looks at the r.h.s. of (A.l) as an algebraic 
equation for p^, so that the equation is a quartic in this variable, then it is possible to show that 
it has always four real roots. This allows to arrange them in order of crescent magnitude and 
then to infer from this fact some properties of the solution of (A.l). The Casimirs and integrals 
are fixed if one fixes the initial conditions, so we can assume that the dynamical variables in 
these quantities are specified by their initial values, say at r = 0. With this in mind, by posing 
P^(t) = a; we rewrite (A.l) as: 

(x^ - 2C2) (x^ - 2Hi) - {Mx -Cif = fix). (A.3) 

If we can find five distinct points where f{x) changes its sign, then the l.h.s. of equation (A.3) 
has four real roots. These points are collected in Table 1. 



'^See for an example [3, page 174]. 
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Table 1. The changes in the sign of f{x). 



X 


fix) 


+00 
VC2 + Hi 

pHo) 

-VC2 + H1 

—00 


>0 

<o 
>o 
<o 

>0 



Note that iff the initial condition are such that J^(0) = J^(0) = p^{0) = p^{0) = 0, then 
p^(0) is equal to one of the two points ±^^2 + Hi. But in this case the four roots are x = p^(0), 
X = J^(0) — p^(0), X = — J^(0) — p^(0), with X = p^{0) a double root, so that also in this case 
there are four real roots. So in general the l.h.s. of (A. 3) has four real roots, with at most three 
equals and with at least one negative (for obvious reasons we do not consider the trivial case 
when all the dynamical variables are initially equal to zero). Given the reality of the roots, it is 
possible now to sort them in an increasing order of magnitude so that by labelling with a, b, c, d, 
we can assume a > b > c> d. Note also that p^{0) lies in the interval (c, b). Equation (A.l) can 
be written then as: 



ipr = iP-a)ip'-b){p'-c)ip'-d). 



(A.4) 



The integration of (A.4) is reduced to a standard elliptic integral of the first kind by the sub- 
stitution [13] z"^ = (^ :}^J\\('lZc\ ■ After some algebra we obtain: 

b — fi'^a sn(ti + rjT, k)'^ 



pHt) 



where 



/^ 



1 — /x^ sn(u + 7]T, kY 

2 _ {a-d){b-c) 



(A.5) 



k' 



c){b-d) 



{a-c){b-d)' 



sn 




(b-p^(0))(a 
a-p3(o))(6 



-,k 



The symbol "sn" denotes the Jacobi elliptic function of modulus k. As can be seen by the 
expansion of p'^{t) in the neighborhood of r = 0, the sign of t] has to be chosen according to 
the sign of p3(o), i.e. sgn{r]) = -sgn{p^{0)) = -sgn(j2(0)p^(0) - .P{0)p^{0)). Note that p^{t) 
is bounded in the set {c,a), that is c < p^{t) < a. Having (A.5) it is a simple matter to write 
down the expressions for the other dynamical variables: 



,{t) = v'2C2 - (p3{T))2, a{T) = v'2ffi - (p'(r))2, cos(*(t)) 



/M^/(0) + f(^ 



Ci - Mp^{t) 
s{T)a{T) 



" ^''^^-\p^{z)cos{ij{z)) - ^M 1 dz. 
z) p 



The continuous flow is in general quasi-periodic. One can ask what are the conditions under 
which the flow becomes periodic. Note that p'^(t), s(r), o"(r) and cos('(/'(r)) are all periodic 
functions with the same period (obviously we are interested only in real periods). Such period 
is given by 



T 



2K{k) 



(A.6) 
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and depends only on the initial conditions and not on the parameters a and /3 of the model. 
In (A. 6) K{}i) is the complete elliptic integral of first kind [13], k and r/ are as given in (A. 5). 
If the function /(t) is such that 

f{T + nT) = fiT) + 2rmT, {n,m) ^1? (A.7) 

then it isn't difficult to see that the motion is indeed completely periodic. In this case the 
trajectories of the moving point (p^(T),p^(T),p^(r)), constrained on the two-sphere with constant 
radius given by r^ = p(t) • p(t), close. The condition (A.7) is equivalent to the following 
constraint on the ratio of the parameters a and /3: 

^ MT 
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